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SUMMARY 

New  constitutive  relations  for  earth  materials  are  aval lable  which 
appear  to  have  theoretical  advantages  over  previous  mathematical  models. 

The  present  work  applies  the  new  model,  which  is  adapted  to  a  tonal Ite,  to 
numerical  computations  of  spherical  wave  propagation.  Parameters  such  as 
attenuation  rates,  stress  trajectories  and  displacement  of  the  cavity  wall 
are  compared  with  similar  parameters  obtained  from  numerical  computations 
using  different  mathematical  models  of  the  same  material.  Except  for  the 
displacement  of  the  cavity  wall,  the  new  model  leads  to  results  which  are 
similar  to  those  found  by  the  most  advanced  of  the  previous  models. 
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SECTION  1 
INTRODUCTION 


New  constitutive  equations  are  available  for  Cedar  City 
tonal i te  which  appear  to  have  theoretical  advantage  over  some 
previous  models  (Reference  1).  The  present  analysis  extends  the 
work  of  Reference  2  by  using  this  mode)  to  calculate  spherical 
wave  propagation.  A  comparison  is  made  between  the  models  con¬ 
sidered  In  Reference  2  and  the  new  "capped"  model,  on  the  basis 
of  attenuation  rates,  stress/time  histories,  cavity  growth,  and 
rate  of  energy  absorption.  An  adaptation  of  the  capped  model  is 
reported  which  appears  to  be  stable,  convergent,  and  applicable  to 
both  finite  element  and  finite  difference  methods  in  continuum 
mechanics. 


The  present  calculations  are  performed  in  the  manner 
described  in  Reference  2.  The  finite  element  method  is  adapted  to 
spherical  geometry.  The  mathematical  mode)  of  Cedar  City  tonal ite 
reported  in  Reference  1  is  modified  such  that  the  shear  modulus, 
hydrostat,  and  yield  criteria  are  as  close  to  those  used  in  Ref¬ 
erence  2  as  the  different  forms  of  the  various  models  permit.  The 
tonal ite  is  considered  to  be  an  infinite  medium  surrounding  a 
cavity  which  contains  a  sphere  of  chemical  explosive.  Detonation 
of  the  explosive  is  represented  by  applying  pressure  to  the  sur¬ 
face  of  the  cavity  which  leaves  with  time  in  a  manner  similar  to 
that  measured  by  Physics  International  Company  (Reference  3)  in 
physical  experiments.  The  peak  pressure  at  the  cavity  wall  is 
31.5  kb. 


The  results  of  the  calculations  are  presented  as  stress/ 
time  histories  and  stress/strain  relations  at  various  ranges  and 
as  rates  of  attenuation  of  the  peak  radial  stress.  The  main 
question  to  be  investigated  is  whether  using  the  device  of  strain 
hardening  to  Incorporate  volumetric  hysteresis  dramatically 
affects  the  wave  shape,  the  attenuation  rate,  or  the  cavity 
growth.  The  findings  indicate  that,  at  least  in  spherical  wave 
propagation,  the  "capped"  model  leads  to  approximately  the  same 
results  as  found  using  the  best  of  the  previous  models. 


--  •  ■  -■ 
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SECTION  2 

FINITE  ELEMENT  ADAPTATION  OF  THE  CAPPED 
MODEL  FOR  ROCK 


The  model  consists  of  an  Ideally  plastic  modified 
Drucker-Prager  yield  criterion, 


fi(v  fl)  ■  ° . <2-’> 

and  a  cap. 

f 2 (J i  *  yjj.  <)  "  0 . (2-2) 


where  is  the  first  stress  invariant,  J2  is  the  second 
Invariant  of  the  stress  deviators,  and  k  is  a  given  function  of 
the  stress-strain  history  of  the  rock.  A  variable  bulk  modulus 
model  which  is  a  function  of  J)  and  a  constant  shear  modulus  are 
used.  The  plastic  potential  flow  rule  is  used  to  determine  plas¬ 
tic  strain  rates  on  both  fj  and  f2> 

YIELD  FUNCTIONS— FRACTURE  SURFACE  AND  CAP 

From  experimental  data  on  Cedar  City  Tonalite,  Sandler 
and  Dimaggio  proposed  the  following  functions  for  fj  and  f^: 


0.0029  J, 

yjJ  -  152  lA5e  -  0  .  . 

.  .  .  (2-3) 

(J,  -  Jc)2 

y  -  Q  -  0 . 

.  .  .  (2-4) 

where  Jj  and  J2  are  expressed  in  ksi.  The  quantities  Jc,  R 
and  Q  are  functions  of  the  cap  parameter  k.  With  reference  to 
Figure  2-1,  Jc  is  the  value  of  Jj  at  the  center  C  of  the 
elliptical  cap,  R  is  the  ratio  of  the  major  to  minor  axes,  and 
Q  is  the  value  of  J2  on  the  cap  when  ■  Jc.  To  ensure  tan- 
gency  at  the  junction  F  of  f^  and  f2»  these  quantities  are 
determined  as  follows: 


2 


•  '  • 
.  I 


AlA 


The  c^p  parameter  k  Is  defined  by 
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i  i 


£  •  f,M  +  KP2)2  +  (de33)2  ;  •  -  (2*,0> 


where  de^*,  de^*  and  d?^  are  the  principal  components  of 
the  plastic  strain  Increments. 


.P* 


BULK  AND  SHEAR  MODULI 

i 

The  bulk  modulus  B  Is  a  function  ofime^n  stress  as 


fot lows: 


&3J, 


B  -  B,  (l  -  B2B3  ’) 


(2-11) 


The  shear  modulus  ,  y  is  a  constant.  , 

,  i 

i  i 

1  t  i 

.  i  1 

y  ■  constant  . . ;  .  .  (2-12) 


i  i 


INCREMENTAL  STRESS/STRAIN  RELATIONSHIP  ,• 

*  .  1 
For  small  strain  Increments,  the  Incremental  Stress 
tensor,  do|j,  is'  related  to  total  strain  Increment  dejj  and 
tha  plastic  strain  increment  de??  as  follows: 


*J , 


»  I 

d°IJ  -  * Kk  -  <)slj  +  2u(dEl  J  -  <})  ■  ■  <2-'3) 


where'  X  and  y  are  Lame  Is  constants,1  5|j  is  the  Kronecker 
delta,  and  the  sum  of  the  elastic  and  plastic  parts  Is  the  total 
strain  Increment. 


The  yield  function  can  be  expressed  3s 


f  (0  j  j ,  k)  1  -  0 


(2-14) 


where  <  is  the  cap  parameter  defined  In  Equation!  2-10.  it  is  to 
be  noted  here  that  the  failure  envelope  is  a  function  df  the 


.1 


1  1 


1  > 


1 


>**■?■■»  •  ••  * 


t 

stresses  and,the  cap  parameter.  An  associated  flow  rule  is  used 
in  conjunction  with  the  capped  yield  surface,  i.e., 


dcf.  -  A#- 
ij  i  3o j  j 


(2-15) 


where  A  Is  a  sdalar  to  be  determined. 

During  Inelastic  deformation,  the  total  differential  of 
the  yield  function  Is  equal  to  zero.  Thus, 


df  ■  f ,  .do. .  +  f  die  ■  0 

,  U  IJ  < 


.  (2-16) 


where 


and 


3f 


IJ  3o. 


3f 

"  3k  • 


The  expression  dx  can  be  obtained  from  Equation  2-10  as 


*'  -  -fl  +  (de33 j' . <2"» 


Substituting  Equation  15  in  Equation  2-17  and  Equation  2-13,  the 
following  relations;  are  obtained. 


and 


d<  -  -A  ft  ^7+  f22  +  f33 


do 


ij 


U 


ij  T  “1WIJ, 


(2-18) 


X  dekk6 i j  +  2u  deiJ  "  A(Xfkk6lj  +  2y6ij)  •  ^2"19^ 


1 


I 


1 

.4.1 


•••  -X-Jw*.  «;»**• 


r\+  tj 


f 


I 


where 


X 


(2-22) 


Substitution  of  Equation  2-21  into  Equation  2-19  yields  the  desired 
Incremental  stress/strain  relationship.  Adjusting  the  tensorial 
shear  strain  e|i  to  the  engineering  shear  strain  y | j ,  the 
following  stress/strain  relationship  is  obtained. 


(2-23) 


6 


8'MMrtw  we  w~r  *»  >  - 1 


AjA 

where 
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F  -  Jl_  + -2L  + -ii- 
30rr  3oee  3ozz 


3f 


rz  3o 


rz 


Y  ,  etc.  ■  X  F  +  2y 


3f 


0  - 


30 


rr 


x  f2  ♦  2U(f2rr  *  f2e  ♦  f2J 


+  f 


ifk  V 


frr  ♦  fee  ♦  fL 


Equation  2-22  Is  the  Instantaneous  incremental  stress/ 
strain  relation  for  small  strain  Increments.  An  Iterative  scheme 
(Figure  2-2),  similar  to  the  modified  Euler  formula,  Is  employed 
to  obtain  stress  Increments  for  f'^lte  Increments  of  strain. 


FIGURE  2-2.  ITERATION  SCHEME  TO  EVALUATE  STRESSES 
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The  procedure  to  obtain  the  stress  increment  during  a  step  is 
described  below. 


Mktl  -  r{Iclk  ♦  ■c>k+i}*4cU . .  •  (2-w 


where  [Clk  and  [C]k+i  are  the  instantaneous  stress/strain 
relationships  at  the  beginning  and  the  end  of  the  step,  respec¬ 
tively.  Since  [C]k+1  depends  on  the  stresses  at  the  end  of 
the  time  step,  successive  corrections  to  [C]^  and  to 
{Ao)k+j  are  made  until  the  following  function  of  stress  incre¬ 
ments  reaches  a  stationary  value  within  a  prescribed  tolerance: 


E  -  [AoJ  +  |Ao2|  +  |AOj| 


(2-25) 
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SECTION  3 
NUMERICAL  RESULTS 


Calculations  were  performed  in  which  a  pressure/time 
history  is  applied  to  the  surface  of  a  spherical  cavity  In  an 
infinite  medium.  The  finite  element  grid  that  represents  the 
medium  and  the  cavity  is  shown  In  Figure  3“1.  The  pressure/time 
history  applied  to  the  surface  of  the  cavity  Is  also  shown. 

A  series  of  eight  calculations  has  been  performed  to 
study  the  Influence  on  spherical  wave  propagation  of  such  material 
parameters  as  volumetric  hysteresis  and  plastic  di latency.  The 
purpose  of  this  study  was  to  determine  which  aspects  of  a  material 
model  most  strongly  affect  calculations  of  wave  propagation  and 
thus  require  the  most  complete  experimental  determination.  The 
first  seven  calculations  (Cases  1  through  7)>  which  explore  a 
range  of  mathematical  models  from  linearly  elastic  through  vari¬ 
able  modulus  with  plasticity,  are  reported  In  Reference  2.  The 
eighth  calculation  (Case  8),  reported  below  considers  an  elastic/ 
plastic  model  with  a  strain  hardening  cap  whose  purpose  Is  to 
provide  plastic  volumetric  compaction  within  the  general  frame¬ 
work  of  plasticity  theory.  Table  3~i  summarizes  the  material 
properties  considered  In  Cases  through  8. 


•  FIGURE  3-1.  MESH  SIZE  AND  LOADING  USED  IN  CALCULATIONS 
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TABLE  3-1.  SUMMARY  OF  MATERIAL  PROPERTIES  AND  LOADING  USED  IN 
SPHERICAL  WAVE  CALCULATIONS 


Case 

Description 

1 

Linear  elastic  (no  yield  condition  imposed); 

B  -  Bq  -  1.2  x  106  psl;  G  -  1  x  106  psl 

2 

Variable  modulus  (no  yield  condition  imposed); 

B  defined  In  Table  3-2;  G  -  1  x  10°  psl 

3 

Variable  modulus  as  In  Case  2  with  von  Mtses 
yield  criterion  (a  ■  0,  c  ■  76,125  psl) 

A 

Variable  modulus  as  In  Case  4  with  Coulomb 
yield  criterion  (a  ■  -0.167,  c  ■  1450  psi)  at 
p  >  -14. 5  kb  and  von  Mlses  criterion  (o  ■  0, 
c  ■  76,125  psi)  at  p  <  - 1 4 . 5  kb.  Plastic 

Potential  flow  rule  Is  used. 

5 

Same  as  Case  4,  except  that  Prandl t-Reuss  flow 
rule  is  used. 

6 

Same  as  Case  4,  except  that  impulse  is  increased 
by  a  factor  of  2. 

7 

Same  as  Case  2,  except  that  impulse  is  Increased 
by  a  factor  of  2. 

8 

Yield  criterion,  bulk  and  shear  moduli  similar 
to  Cases  4  and  5i  except  that  capped  model  is 
used  to  represent  volumetric  hysteresis. 

NOTE:  Cases  1  through  7  are  previous  analyses;  Case  8  Is  the 
present  analysis.  / 


The  yield  criteria  and  hydrostatic  stress/straln  prop¬ 
erties  of  the  material  model  in  Cases  4,  5,  and  8  are  similar  as 
indicated  in  Table  3-2  and  Figures  3-2  and  3~3>  However,  the 
mathematical  formulation  of  the  capped  model  (Case  8)  differs 
significantly  from  that  of  the  other  two.  The  models  In  Cases  4, 
5,  and  8  and  known  to  be  similar  only  for  certain  stress  paths. 
Differences  may  occur  when  the  material  is  subjected  to  other 
stress  paths.  The  subject  which  is  investigated  below  Is  whether 
the  differences  in  mathematical  formulation  lead  to  differences 
in  the  response  to  spherical  wave  propagation. 
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FIGURE  3-3.  BULK  MODULUS  AND  HYDROSTAT  USED  IN  CASE  8  ARE  PRESCRIBED 
BY  DIFFERENT  FORMULAS  BUT  HAVE  NEARLY  THE  SAME  VALUES  AS 
IN  PREVIOUS  CASES 


(b)  HYDAOSTAT,  CASES  k,  5  AMD  8 


FIGURE  3-3.  (CONTINUED) 
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SECTION  4 

RESULTS  OF  CALCULATIONS 


STRESS/TIME  HISTORIES 

The  radial  and  circumferential  stress/time  histories  for 
Cases  4,  Si  and  8  are  shown  in  Figure  4-1.  Compressive  stresses 
are  defined  as  negative. 

Cases  4  and  8  are  similar  in  that  the  circumferential 
stresses  remain  compressive  throughout  the  duration  of  the  calcu¬ 
lation.  They  differ  from  Case  5  in  this  respect  for  the  following 
reason.  The  yield  criteria  have  a  Coulomb  portion  which  allows 
plastic  deformation  to  occur  at  all  ranges  considered  in  the  three 
cases.  The  materials  of  Cases  4  and  8  exhibit  dl latency,  or 
volumetric  expansion,  which  tends  to  induce  lateral  compressive 
stresses.  This  tendency  is  partly  counteracted  by  outward  radial 
displacement,  which  contributes  to  lateral  tensile  stresses. 

Oi latency  dominates  in  Cases  4  and  8,  however,  in  that  tensile 
hoop  stresses  do  not  develop.  In  Case  5,  no  di  latency  is  allowed, 
and  hence  tensile  stresses  develop. 


STRESS  TRAJECTORIES 

The  stress  history  in  terms  of  or  versus  oQ,  o  is 
referred  to  below  as  a  stress  trajectory.  Trajectories  at  ^ 
r/r  *  1,  2,  3,  and  4  are  shown  in  Figures  4-2  through  4-5.  The 
loading  is  initially  elastic  up  to  about  ■  -7  Kb.  Loading 
continues  on  the  cap  up  to  the  peak  values  of  or  and  erg,  ox. 
Unloading  is  initially  elastic  in  the  region  where  oscillations 
occur.  As  unloading  continues,  the  stress  point  again  encounters 
the  cap.  At  r/rQ  ■  3  and  4,  the  stress  point  remains  on  the  cap 
throughout  the  duration  unloading.  At  r/r0  »  1,  the  stress  point 
spends  only  a  brief  time  on  the  cap  and  then  moves  onto  the 
fracture  surface.  At  r/rQ  ■  2,  the  stress  point  is  primarily  in 
the  elastic  region  and  on  the  cap  and  moves  onto  the  fracture 
surface  only  near  the  end  of  the  calculation. 

During  initial  loading,  the  cap  moves  rapidly  away  from 
the  origin.  In  this  phase  there  is  at  some  stages  plastic  volu¬ 
metric  compaction  and  at  other  stages  plastic  volumetric  expansion. 
There  is  also  some  plastic  shear  strain.  During  unloading  on  the 
cap,  the  sense  of  the  plastic  volume  strain  is  exclusively  expan¬ 
sion.  The  motion  of  the  cap  during  unloading  is  small  compared 
with  that  during  loading.  This  is  because  changes  in  the  cap 
parameter  depend  on  fj  (see  Equation  2-10),  which  tends  toward 
zero  as  the  stress  point  approaches  the  fracture  surface. 
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DILATENCY  AND  t/y 

In  Reference  2,  an  analysis  was  made  of  Inelastic  energy 
absorption  at  several  ranges.  In  Case  4  it  was  found  that  di la¬ 
tency,  which  arises  from  using  the  plastic  potential  flow  rule 
combined  with  yield  criterion,  produces  a  net  volumetric  expansion. 
It  is  suggested  that  one  effect  of  di  latency  is  to  reduce  the 
unloading  wave  velocity,  and  that  this  is  why  the  peak  radial 
stress  attenuates  more  slowly  in  Case  4  than  in  Cases  2  and  5  (no 
di  latency). 

To  investigate  this  subject  for  the  capped  model,  plots 
of  pressure  versus  excess  compression  (P/p)  and  shear  stress 
versus  shear  strain  (t/y)  were  prepared  and  are  shown  in  Fig¬ 
ures  4-6  and  4-7.  A  main  finding  of  Case  8  is  that  there  is  net 
di latency  at  r/rQ  ■  2  and  3  (Figure  4-6(c)),  just  as  there  is 
in  Case  4.  Although  It  is  difficult  to  define  the  term  "amount  of 
di latency,"  there  would  appear  to  be  less  in  Case  8  than  in  Case  4. 
If  this  is  true,  of  the  three  cases  in  Figure  4-6,  Case  4  has  the 
greatest  amount  of  di  latency,  Case  8  is  intermediate,  and  Case  5 
has  none.  This  is  a  useful  observation  in  interpreting  the 
attenuation  rates. 


ATTENUATION  RATES 

The  rates  at  which  peak  radial  stresses  attenuate  with 
range  in  Cases  4,  5,  and  8  are  shown  in  Figure  4-8.  Beyond 
r/rQ  -  2.5  the  attenuation  rate  in  Case  8  is  intermediate  between 
Cases  4  and  5.  Thus,  as  in  the  previous  cases,  the  attenuation 
rate  appears  to  be  correlated  with  the  amount  of  di  latency  and 
hence  the  speed  of  unloading  signals. 


GROWTH  OF  CAVITY 

The  displacements  of  the  cavity  wall  for  Cases  4,  5, 
and  8  are  shown  in  Figure  4-9.  The  relative  magnitudes  of  the 
displacements  reflect  important  aspects  of  material  behavior. 

Up  until  5  x  10”°  sec,  the  displacement  in  Cases  4  and  5  are 
similar.  However,  as  soon  as  unloading  proceeds  to  where  the 
stress  point  is  on  the  yield  surface,  di latency  develops  In 
Case  4  and  does  not  develop  in  Case  5.  Since  di  latency  tends  to 
retard  cavity  growth,  the  cavity  grows  more  rapidly  in  Case  5 
than  in  Case  4. 

The  reason  for  greater  cavity  growth  in  Case  8  is  partly 
due  to  the  fact  that  plastic  shear  strains  occur  on  loading  as  a 
result  of  using  the  plastic  potential  flow  rule  in  conjunction 
with  the  cap.  This  cannot  be  the  full  explanation,  however.  The 
reason  cited  above  should  affect  displacements  more  at  early  than 
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FIGURE  4-7.  SHEAR  STRESS/STRAIN  CURVES  AT  r/r  -  2  AND  r/r  -  3 
FOR  CASES  4,  5  AND  8  °  ° 
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FIGURE  4-8.  ATTENUATION  OF  PEAK  FIGURE  4-9.  DISPLACEMENT  OF  CAVITY 
RADIAL  STRESS,  WALL,  CASES  4,  5  AND  8 

CASES  4,  5  AND  8 


at  late  times.  Figure  4-9  shows  that  displacements  in  Cases  5  and 
and  8  are  identical  at  early  times  and  do  not  begin  to  differ 
until  after  20  usee. 


EFFECTIVENESS  OF  THE  ITERATIVE  TECHNIQUE 

The  reasons  for  resorting  to  the  iterative  technique 
described  in  Section  2  is  that  stress/strain  coefficients  and 
the  cap  parameters  are  implicit  functions.  Since  there  are  many 
possible  iterative  techniques  and  the  one  presently  used  is  some¬ 
what  complicated  and  time-consuming,  it  should  at  least  be  demon¬ 
strated  to  be  accurate.  Two  demonstrations  are  made  below.  First, 
a  sequence  of  shock  loading  stress  states  is  examined  at  the  end 
of  each  time  step  to  see  whether  they  lie  on  the  appropriate  cap.* 
Second,  the  time  history  of  the  cap  parameter  for  the  element 
adjacent  to  the  cavity  is  studied  to  determine  whether  it  varies 
with  the  integration  time  step  (it). 


*This  is  a  more  stringent  test  than  showing 
on  the  fracture  surface  (which  can  also  be 
cap  moves,  whereas  the  fracture  surface  is 


the  stresses  to  lie 
done)  because  the 
constant. 
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In  Case  8(a)  (it  ■  5  x  10"8  sec),  the  value  of  the  cap 
parameter  <  in  element  t  (adjacent  to  cavity)  increases  from 
zero  to  about  330  psi.  The  caps  corresponding  to  shock  loading 
and  early  unloading  are  defined  by  0  <  k  <  175  psi  and  are 
shown  in  Figure  4-10.  The  number  of  time  steps  between  successive 
caps  varies  from  1  (0  <  <  <  86  psi)  5  to  7  (145  psi  <  k  s  175  psi). 
The  maximum  change  in  k  in  any  one  step  is  about  120  percent 
(k  ■  31  psi  -*•  k  -  70  psi).  Figure  4-10  shows  that  the  stress 
state  lies  exactly  on  the  cap  up  until  k  -  137  psi  and  begins  to 
drift  off  at  higher  values.  The  reason  for  this  drifting  Is  sug¬ 
gested  in  Figure  4-10,  where,  for  k  >  125  psi,  part  of  the  cap 
lies  above  the  fracture  surface.  The  model  should  not  behave  in 
this  way.  The  flaw  in  the  model  appears  to  lie  in  the  functions 
contained  in  Equation  2-6.  While  R*  increases  very  rapidly  with 
Jjp,  the  term 


NOTE: 

NUMBERS  REFER  TO  VALUE  OF  CAP  PARAMETER  <  IN  PSI.  ^/DENOTES 
THE  APPROXIMATE  POSITION  WHICH  THE  STRESS  POINTS  SHOULD  HAVE 
ON  ITS  CORRESPONDING  CAP  IF  ITERATION  TECHNIQUE  WERE  EXACT. 


J,,  KSI 


FIGURE  4-10.  RELATIVE  POSITIONS  OF  SUCCESSIVE  STRESS  POINTS 
AND  CORRESPONDING  CAPS,  ELEMENT  1,  CASE  8-a 
(At  -  5  x  10"8  SEC) 
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should  tend  toward  zero  at  the  same  rate.  Thus  /Q  (see  Figure  2-1) 
would  tend  toward  the  von  Mises  limit  of  the  fracture  surface. 
Unfortunately,  this  balance  is  not  achieved  and  /Q  exceeds  the 
von  Mises  limit.  The  stress  point  eventually  becomes  constrained 
by  the  fracture  surface  and  thus  drifts  away  from  the  cap.  The 
cap  becomes  larger  and  larger  without  any  control  being  imposed 
by  the  fracture  surface. 

In  spite  of  the  apparent  success  of  the  present  itera¬ 
tive  technique,  some  difficulties  remain  as  illustrated  in  Fig¬ 
ure  4-11.  In  this  figure  the  time  histories  of  ic  in  Element  1 
are  shown  for  Cases  8(a)  and  8(b).  The  only  difference  between 
the  calculations  is  the  integration  time  step  which  is  5  x  10~8  sec 
in  8(a)  and  1.5  x  i10"7  sec  in  8(b).  Case  8(a)  uses  the  same  time 
step  as  the  previous  seven  cases  and  is  considered  to  be  superior 
to  Case  8(b).  One  of  the  reasons  Case  8(a)  is  considered  to  be 
superior  is  that  the  cap  parameter  is  not  required  to  undergo 
such  large  changes  in  one  time  step.  Figure  4-11  shows  that  the 
present  iterative  procedure  works  about  equally  well  for  the 
large  as  for  the  small  time  step  in  the  early  times.  As  the  cap 
approaches  the  yield  surface  however,  and  the  model  develops  the 
flaw  discussed  above,  the  values  of  k  tend  to  diverge.  This 
has  some  effect  on  response  parameters  of  interest  as  reflected 
In  the  different  displacements  of  the  cavity  wall,  Figure  4-12. 


FIGURE  4-11.  CHANGE  IN  SURFACE 
(r/r0  -  1.015)  CAP 
PARAMETER  k  WITH 
TIME  AT  CAVITY 


FIGURE  4-12.  DISPLACEMENT  OF 

CAVITY  WALL  IN  CASE  8 
FOR  TWO  DIFFERENT 
TIME  STEPS 
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SECTION  5 

SUMMARY  AND  CONCLUSIONS 


Under  conditions  of  spherical  wave  propagation,  the 
capped  model  of  rock  is  similar  to  the  elastic,  perfectly  plastic 
models  of  Reference  2  with  respect  to  wave  shape  and  attenuation 
rate  of  peak  stress.  The  amounts  of  energy  absorbed  in  shear 
deformation  and  gained  in  volumetric  dilatation  are  also  close  to 
those  found  for  the  models  of  Reference  2.  The  displacements  of 
the  cavity  wall  are  significantly  greater  in  the  capped  model  than 
those  in  the  previous  models.  The  reason  for  the  difference  is 
not  fully  explained.  Potential  numerical  difficulties  associated 
with  keeping  the  stress  point  on  the  cap  or  on  the  fracture  sur¬ 
face  during  plastic  deformation  are  noted.  These  are  overcome  in 
the  present  work  by  means  of  an  iterative  technique  which  uses  an 
Euler  formula  for  predicting  and  correcting  the  inelastic  stress/ 
strain  coefficients  until  a  stationary  value  of  an  arbitrary 
stress  function  is  reached.  A  defect  in  the  operation  of  the 
present  capped  model  is  noted,  whereby  it  is  possible  for  the  cap 
to  get  outside  the  fracture  surface. 

The  main  conclusions  of  the  present  study  are  as  follows: 

a.  In  applications  to  spherical  wave  propagation,  the 
capped  model  gives  results  which,  in  most  respects, 
agree  with  the  results  from  previous  models. 

b.  An  iterative  technique  which  uses  an  Euler  formula 
to  predict  and  correct  tangent  moduli  appears  to 
be  well-suited  to  finite  element  applications 
when  the  model  contains  implicit  functions  of  the 
current  plastic  strain  increments  and  the  total 
stresses.  The  method  can  also  be  used  for  finite 
difference  applications. 
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